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Most Supermassive Black Holes must be Rapidly Rotating 

M. Elvis 1 2 , G. Risaliti 1,2 & G. Zamorani 3 

ABSTRACT 

We use the integrated spectrum of the X-ray background and quasars Spectral 
Energy Distribution to derive the contribution of quasars to the energy output 
of the Universe. We find a lower limit for the energy from accretion onto black 
holes of 6%, of the total luminosity of the Universe and probably more, with 
15% quite possible. Comparing these values with the masses of black holes in 
the center of nearby galaxies we show that the accretion process must be on 
average very efficient: at least 15% of the accreted mass must be transformed 
into radiated energy. This further implies that most supermassive black holes 
are rapidly rotating. 

Subject headings: galaxies: active — X-rays: diffuse background 


1. Introduction 

The recent realization that most of quasars are hidden behind large amounts of dust 
and gas (Fabian & Iwasawa 1999, hereafter FI99, Maiolino & Rieke 1995) has increased the 
estimate of the total quasar luminosity of the universe substantially, to a point where it 
is no longer negligible compared with that of stars, and to a level where it fully explains 
the diffuse background seen in X-rays. Since Hubble observations have determined the mean 
mass density of black holes, the ratio of luminosity to mass determines a maximum efficiency 
required. If this substantially exceeds 5% then supermassive black holes must, on average, 
be rapidly spinning. 

FI99 demonstrated that assuming 10% efficiency the X-ray background and the black 
hole mass density were consistent. Here we revisit their argument, carefully examining all 
the assumptions, to produce a strong lower limit on the supermassive black hole efficiency. 
We find that this efficiency must be greater than 15%. Most supermassive black holes must 
then be spinning. 
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2. The total luminosity from quasars 


The hard X-ray background spectrum (XRB, E> 2 keV) is accurately reproduced by 
the integrated emission of AGNs, most of which have to be obscured, (Setti & Woltjer 1989, 
Comastri et al. 1995, Gilli et al. 2001). These synthesis models successfully predict the 
Chandra and XMM-Newton hard X-ray source counts (Hasinger et al. 2001). With this 
increasingly well supported model the X-ray background spectrum can be used, together 
with a quasar SED, to estimate the total energy output of AGNs in the Universe, as shown 
by FI99. 

A key step in the argument of FI99 is the bolometric correction assumed to derive the 
total quasar luminosity from the X-ray intensity. The SED assumed by FI99 is that of Elvis 
et al. (1994, hereafter E94)), obtained from a sample of quasars for which an Einstein IPG 
spectrum was available. Here 1) we derive a different bolometric correction by taking into 
account the statistical bias on the average X-ray to optical ratio derived by the analysis of an 
X-ray selected sample and by excluding the infrared region (which would otherwise double 
count the UV and soft X-ray photons); 2) we discuss the possibility that this ratio decreases 
with redshift; and 3) we complete the Elvis et al. SED with the recent UV SED of Cheng et 
al. 1999. All together, these changes increase by a factor > 2 the estimate of the quasars’ 
total energy output with respect to the value derived by FI99. 

We can estimate the total energy released by quasars with the following two steps: 
(1) we assume that the X-ray background at 30 keV is due to unobscured emission from 
AGNs and then extrapolate the value of the unobscured background down to 2 keV, ^ 2 kev = 
E x I(E)( 2keV), assuming an average X-ray spectrum for quasars; (2) we compute the total 
energy output from quasars,using an average quasar SED. The first step involves a 
number of possible corrections, each of which of the order of 30%, with respect to the simple 
case of a single power law slope adopted by FI99; the second step leads to a bolometric 
correction of a factor 1.6-3 higher than that of FI99. 

1) The spectral complexity of quasars in the X-rays leads to four small uncertainties in 
estimating the unobscured emitted 2 keV X-ray background: 

(a) Photon index-, the value T = 2 used by FI99 comes from the subsamples of low redshift 
optically selected PG quasars (Schmidt & Green 1983, George et al. 2000, Mineo et al. 2000). 
However, Reeves et al. (2000) find T = 1.89 ± 0.05 for all radio quiet quasars observed by 
ASCA (including several high redshift objects), which would decrease T 2 keV by 25%. and 
may be due to a flattening of the average X-ray spectrum of quasars at high redshifts. We 
note though that since this sample is X-ray selected there could be a bias towards hard X-ray 
loud objects. 

(b) Reflection component a reflection component in quasars spectra could decrease the 
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prediction of T 2 keV- However, the effect should be small (< 5%), especially since the reflection 
component is weaker in high luminosity quasars (Mineo et al. 2000). 

(c) High energy cutoff. BeppoSAX observations (Guainazzi et ah 1999, Perola et ah 1999, 
Perola et ah 2000, Nicastro et ah 2000, Risaliti et ah 2001) find exponential cutoffs in AGN 
spectra at energies E> 100 keV. Assuming an average redshift of 2 for quasars making most 
of the XRB, and a cut off energy E=300 keV, the normalization of the continuum powerlaw 
(and so of T 2 keV ) is increased by 35%. 

(d) Extremely High N H objects: quasars absorbed by column densities N H > 10 24 cm -2 
can be significantly absorbed even at a rest frame energy corresponding to 30 keV in the 
observer frame. Moreover, the X-ray emission is completely absorbed if N h > 10 25 cm -2 . 
According to Risaliti et ah (1999) objects with N# > 10 25 cm -2 account for at least 20% of 
low luminosity AGNs at low redshift, and could be more than 30%, while another 10-20% 
have 10 24 cm -2 < N# < 10 25 cm -2 . Therefore, if these percentages remain the same at all 
luminosities, we must increase T 2 keV by at least 20%. 

To estimate ^ 2 kev we used the code of Gilli et ah (2001), which reproduces the XRB 
spectrum and source counts (Akiyama et ah 1999, Fiore et ah 1999, Hasinger et ah 2001). 
This code uses the soft X-ray luminosity function of AGNs and its evolution with redshift 
(Miyaji, Hasinger & Schmidt 2000), and assumes an average X-ray SED for type 1 (unob¬ 
scured) and type 2 (obscured) AGNs. To estimate the total intrinsic luminosity emitted by 
AGNs, we assumed an unabsorbed SED for all AGNs with T = 1.9, Ec=300 keV. These 
parameters predict a relatively low XRB intensity at 2 keV: a much higher value is quite 
possible, given the parameter space allowed by the observational results described above. 
The result we obtain is T 2 keV = 40 keV s -1 cm -2 sr^ 1 (only 5% larger than 38 keV s^ 1 
cm -2 sr -1 , of FI99). Assuming a fraction of completely obscured objects / a bs=0.2 we obtain 
T 2 keV — 48 keV s _1 cm -2 srW The Gilli et ah model predicts an increase of the fraction 
of absorbed objects with redshift, so this value, already a lower limit for low redshift AGNs, 
is even more conservative at high redshifts. 

2) We now use T 2 keV, together with a bolometric correction, to estimate the total 
integrated emission of quasars, At wavelengths from the radio to the UV we can assume 
that the SED by E94 is representative for the average emission of unobscured quasars. 
However, the X-ray emission of quasars in the E94 sample included only those quasars that 
gave a large enough number of photons in the Einstein IPG to construct a well-constrained 
spectrum. Since the exposure times for the optically selected quasars were all similar, the 
E94 sample is biased toward X-ray brighter quasars. This is a sizeable effect. The commonly 
used indicator for the X-ray to optical ratio in quasars is the a (ox) index, defined as a (ox) = 
— log(/^// o )/2.605, where f Q and f x are the monochromatic fluxes at 2500A and 2 keV. 
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Given a sample of AGNs with an observed normalized distribution D(a(ox)) of the 
a (ox) parameter, we can define an effective a (ox) (see Avni et al. (1980) for details) as: 

< fo/fx >= 10 2 - 605x “K)e// = I D(a(ox)) x io 2m5xa ^da(ox) (1) 

Since our starting point is an integrated X-ray emission (i.e. T(A")), the corresponding 
< fo/fx > to be used for the estimate of the appropriate X-ray to optical conversion is 
obtained from Eq. 1 with the parameters typical of X-ray selected samples. Assuming, 
consistently with the data, a Gaussian distribution of a (ox), from Eq. 1 we obtain 

a(ox) e ff =< a(ox) > +3<r 2 (2) 

where < a (ox) > and a 2 are the mean and variance of D(a(ox)). A reasonable value for a 
is ~ 0.16, in agreement with Yuan et al. 1998. Using < a(ox) >= 1.35 (from E94, and in 
agreement with the results from the ROSAT deep surveys (Lehmann et al. 2001), we obtain 
a(ox) e ff = 1.43. 

Reassuringly, we obtain the same result starting from the optically selected sample of PG 
quasars (Schmidt & Green 1983). In this case we have < a(ox) >= 1.55 for the subsample of 
low redshift PG quasars (Laor et al. 1997). If a(ox) has a gaussian distribution, the relation 
between the observed values of a (ox) in flux limited X-ray and optically selected samples 
(assuming the same parent distribution) is 

< a (ox) > op t ~ < a (ox) > x = 6 cr 2 (/3 — 1) (3) 

where (3 is the differential slope of the counts as a function of flux (see Francis et al. (1993) 
for a similar discussion in a different context). Assuming the Euclidean slope f3 = 2.5, the 
observed < a (ox) > for this optically selected sample is approximately equivalent to our 
adopted < a (ox) > for an X-ray selected sample and therefore consistent with a(ox) e ff = 
1.43. The use of this value increases the bolometric correction by a factor 1.6, with respect 
to the value a (ox) =1.35 used by FI99. 

Most of the XRB originates at substantial redshifts, z~ 1 — 2. A study of the optically 
selected quasars observed by ROSAT (Yuan et al. 1998) suggests a possible increase of 
a(ox) with redshift or luminosity to an average a (ox) =1.68 at z=2. If this increase is due 
to a luminosity effect, since most of the XRB is made by objects of moderate luminosity 
(L(2-10 keV)< 10 45 erg s -1 ), the safest assumption is to assume the < a(ox) > value 
obtained with the complete low redshift sample, which is free from biases and samples well 
the relevant luminosity range. If, instead, the increase in < a(ox) > is mainly correlated 
with redshift, then it would affect the bolometric correction. We can estimate an upper limit 
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to this correction by using the value < a(ox) >= 1.68 at z=2 (Yuan et al. 1998). The 
corresponding < a(ox) > e ff would be ~ 1.55 and the bolometric correction would be 3.3 
times higher than that estimated by FI99. 

3) Final source of uncertainty is the EUV SED of quasars. The SED of E94 does not 
cover the energy range between the H edge (E=13.6 eV) and E~ 0.15 keV. Here we use 
the UV SED obtained by Zheng et al. (1999) from a large sample of radio quiet quasars 
observed with HST in the 350A-912A range. From 350A to 0.15 keV we assume a power law 
spectrum (the details of the SED in this band affect the final bolometric correction by less 
than 5%). The choice of this UV SED is, once again, conservative in the sense to predict 
the lowest possible bolometric correction. For comparison, the SED estimated by Mathews 
& Ferland (1987), based on the equivalent width of high ionization emission lines, leads to a 
~ 25% higher bolometric correction. 

Using the lower limit at 2 keV computed above and assuming a(ox) e ff to be in the 
range of 1.43-1.55 we obtain the integrated emission of AGNs: \I/t > (3.6 — 8) nW m 2 sr 1 , 
where the first value is likely to be a lower limit. This estimate has to be compared with the 
value of ~ 2 nW m _2 sr _1 , obtained by FI99 (actually FI99 obtain a value of 3 in the same 
units, but, as they say in the Discussion, this value has to be lowered by 1/3, in order not to 
count the UV emission twice, since the IR radiation is due to reprocessing of UV emission). 

Converted into an energy density this gives: 

ZLyrtf/ 

Ut =-~ (1.5 — 3.4) x 1CT 15 erg cm -3 (4) 

c 

If we assume that this energy comes from accretion with efficiency e, the mass density in 
black holes is: 

p. = ^ ~ (7.5 - 16.8) x 10 4 e _1 M 0 Mpc” 3 (5) 

e c 2 

where we assume that the peak of AGN emission is at z=2, in agreement with the XRB 
synthesis models. 

Our results can be compared with estimates of the total luminosity in the Universe from 
both stars and supermassive central black holes in galaxies. 

The total luminosity of the Universe has been estimated from the diffuse background 
from sub-mm to UV wavelengths. Madau & Pozzetti (2000) estimate that the diffuse back¬ 
ground between 0.1 and 7 pm has an integrated intensity ^uo = 17 ± 3 nW m” 2 sr” 1 . 
The bolometric emission in the far infrared, H/jr, between 1000pm and 7pm, derived from 
the IR background, is ~ 40 nW nr” 2 sr” 1 (Franceschini et al. 2001). The stellar produced 
background outside these spectral regions is negligible. Therefore, the total energy in the 
background, due to stars and quasars combined, is T T ot ~ 57 nW m” 2 sr” 1 . This means 
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that quasars contribute at least 6% to the total luminosity of the Universe. This is a lower 
limit only dependent on onr knowledge of the SED of quasars. Since we obtained our es¬ 
timate through a chain of conservative assumptions, a significantly higher contribution is 
highly probable. A value ~ 15% could apply if the redshift evolution of the optical to X-ray 
ratio is real, and up to 20-25% if we release our conservative assumptions on quasar spectra 
and the fraction of Compton thick objects (for example, a photon index T = 2 instead of the 
adopted T = 1.9, and a fraction of quasars with N# > 10 25 cm -2 f=30% instead of f=20% 
are quite possible, if not favored by the hard X-ray data now available). 


3. Accretion efficiency in AGNs 

The mass density in central black holes has been estimated by Salucci et al. (1999) 
to be ~ 4 x 10 5 M q Mpc -3 using a quasar evolution model and starting from X-ray 
counts. These results are in agreement with a recent work (Merritt & Ferrarese 2001) on 
the relation between central black hole mass and the velocity dispersion of the bulge (best 
estimate p, ~5x 10 5 M 0 Mpc -3 ) 

Combining these estimates with Equation 1 gives a value for the average accretion 
efficiency, e. We find e > 0.15. 

Averaged over the history of the Universe our result implies that every black hole is 
radiating at > 15% efficiency. This indicates that the accretion history of the Universe 
cannot be on average silent, and rules out the possibility that a significant fraction of the 
black hole mass accretion has occurred in the ADAF regime. Moreover, this value is well 
above the maximum efficiency of non-rotating black holes, (e^R = 5.6%). Therefore, most 
supermassive black holes must be rapidly spinning. 


This work was supported in part by NASA grant NAG5-4808. 
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